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Non-perturbative effects of constant magnetic fields in a Higgs- Yukawa gauge model are studied 
using the extremum equations of the effective action for composite operators, ft is found that 
the magnetic field induces a Higgs condensate, a fermion-antifermion condensate, and a fermion 
dynamical mass, hence breaking the discrete chiral symmetry of the theory. The results imply that 
for a non-simple group extension of the present model, the external magnetic field would either 
induce or reinforce gauge symmetry breaking. Possible cosmological applications of these results in 
the electroweak phase transition are suggested. 
Symmetry behavior in quantum field theories under the influence of external fields has long been a topic of intensive 
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^_i ' study in theoretical physics Q . In the present paper we are interested in particular in non-perturbative effects produced 
iq , by external magnetic fields in gauge theories with scalars. Our main claim is that for theories with non-simple gauge 
£*\j ■ group, scalar-scalar and scalar-fermion interactions, the magnetic field reinforces gauge symmetry breaking. 

The observation of large-scale galactic magnetic fields in a number of galaxies, in galactic halos, and in clusters of 
OQ ■ galaxies B has recently stimulated a large number of works trying to explain the physical mechanism responsible for 
the origin of these fields. Many of the proposed generating mechanisms have compelling arguments in favor of the 
existence of strong primordial magnetic fields (for a review of cosmological generating mechanisms see H and references 
therein) . Since primordial magnetic fields could play a significant role in particle cosmology, the investigations on the 
theme have recently boomed. In this context, the implications of a magnetic-field-driven gauge symmetry breaking 
— ■ mechanism may be important. 

OS \ Several years before this renewed interest in cosmological magnetic fields, Amb0jrn and Olesen W considered the 

OS ■ electroweak model in the presence of a constant magnetic field. Assuming certain special values of the couplings, 
they obtained a W- and Z-condensate solution forming a lattice of abelian vortex lines for a range of magnetic fields 
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lying between — ^ and mw cos — . At even larger values of the magnetic fields they found that the phase transition 
Qh! to a symmetric phase can be reached at temperatures lower than the critical one at zero field. This result realizes, 
although due to a totally different reason, an old suggestion 0] that large magnetic fields could induce the transition 
from the broken to the unbroken phase in the electroweak system. 

More recently, the ground state of the electroweak theory in the presence of a hypermagnetic field has been investi- 
gated using either numerical or perturbative calculations ||- ||10f . The main motivation of these papers was to study 
the possibility that a hypermagnetic field could allow the realization of baryogenesis within the Standard Model M . 
Even though the original results || for the upper bound of the Higgs mass needed to have baryogenesis in the SM 
were quite optimist, it was quickly realized that higher loop effects pi, [ ]l(| and numerical non-perturbative calcu- 
lations [p| would significantly weaken the transition. Moreover, posterior studies on which certain subtleties of the 
theory- like the magnetic dipole moment of the sphaleron |S[| or ring diagrams contributions to the high-temperature 
effective potential [Q- were taken into account, concluded that albeit the hypermagnetic field strengthen the first 
order character of the phase transition, it is not enough to satisfy the SM baryogenesis condition |flj 

^)>1, (1) 

T c ~ ' K } 

with v (T c ) the Higgs vacuum expectation value (vev) at the critical temperature T c of the electroweak phase transition. 
When a non-perturbative analytic approach is used to study field theories in external magnetic fields, new non- 
trivial effects are found. An important example of these non-perturbative effects is the formation of a chiral symmetry 
breaking fermion condensate < tptp > and of a dynamically generated fermion mass in the presence of an external 
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magnetic field, known in the literature as magnetic catalysis |12|]. This phenomenon, which has proven to be rather 
universal and model independent, has recently attracted a lot of attention |13]- |lq |. 

On normal circumstances massless fermions can condensate and acquire a dynamical mass, but the condensate 
appears only for sufficiently strong coupling between fermions. The new feature when a magnetic field is present is 
that it favors (catalyzes) the symmetry breaking by reducing to the weakest attractive coupling the strength of the 
interaction needed to break the symmetry. The essence of this effect is that the fermions in the lowest Landau level 
(LLL) constitute the effective fermionic degrees of freedom whose dynamics dominates the long wavelength behavior 
of the system. The phenomenon is driven by the fact that massless fermions acquire an energy gap in the presence of 
a magnetic field, but there is no energy gap between the vacuum and the LLL fermions. Then, in the infrared region, 
the dynamics of the LLL fermions dominates the fermion propagator, making it essentially D-2 dimensional. This 
effective dimensional reduction strengthen the fermion pairing dynamics ]l2fl , ]l9| ] giving rise to a fermion condensate. 

It is worth to mention that the phenomenon of magnetic catalysis is not only interesting from a purely fundamental 
point of view, but it has potential application in condensed matter J2Q]- pjj and cosmology [O]. For instance, it 
has been recently speculated that the generation of mass through magnetic catalysis in lower dimensional models 
|pi| , |22| , or in four-dimensional models with boundaries |23| , could be behind the physical mechanism explaining the 
observed scaling of the thermal conductivity in superconducting cuprates with an externally applied magnetic field 
p4f . On the other hand, the magnetic catalysis could influence the character of the electroweak phase transition as 
suggested by the results of ref . |ll| . 

In the present paper we consider a simple model field theory with the aim of investigating in a self-consistent 
way how scalar-scalar and fermion-scalar interactions in the presence of an external magnetic field can influence the 
stability of the vacuum. It is not intended as a realistic theory, but rather as an example of a large class of theories 
with scalar fields, on which dynamical symmetry breaking (either chiral or gauge) can be catalyzed by an external 
magnetic field. In this sense, it could be useful for condensed matter, as well as for cosmological applications. If 
this toy model is extended to include a non-simple gauge group theory, as for instance the electroweak model, the 
results of this paper could provide a scenario on which, in contrast to the effect found by Amb0jrn and Olesen M, 
an external magnetic field could either induce gauge symmetry breaking or reinforce it, in the case it already exists, 
through non-perturbative effects. 

Let us consider the following theory of gauge, fermionic and real scalar fields described by the Higgs- Yukawa 
Lagrangian density 

L = -\f^F^ v + tyfW + g&f^A? - \d^d»y - ±<p* - ^V - \ yV W (2) 



This theory has a U(l) gauge symmetry, 



a fermion number global symmetry 



and a discrete chiral symmetry 



A,, -> A^ + -d^x) 
g 

i, -» e M(x V, (3) 



V -> e'% (4) 



?/>^7 5 V>, v^-Vn^ <p->-<p (5) 



Note that a fermion mass term mipip is forbidden, since it is invariant under (g) and (Q), but not under the discrete 
chiral symmetry (|3|). 

To study the vacuum solutions that could arise in the theory (0) under the influence of an external constant magnetic 
field B, we need to solve the extremum equations of the effective action F for composite operators (2J|, [£6| 

"^-O, (6) 
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Scp c 
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where G(x,x) = <r(x) = (0 | ip(x)ip(x) | 0) is a composite fermion-antifermion field, and <p c represents the vev of the 
Higgs field. Thanks to the discrete chiral symmetry (H), it is enough to consider only one composite field. We choose 



the composite field G(x,x), ignoring the second possible one, n(x) = (0 | ip{x)i^tp{x) | 0) , since the effective action 
can be a function only of the chirally invariant combination p 2 = a 2 + ir 2 . 

The loop expansion of the effective action T for composite operators [£5| , [£6| can be expressed as 

r (G, <p c ) = S (<p e ) - iTrlnW 1 + i-Tr lnD" 1 + i-Tr In A" 1 - iTr [G^ 1 {<p c )G] + T 2 (G, <p e ) + C (8) 

In Eq. (H) G is a constant and S (<p c ) is the classical action evaluated in the scalar vev (Higgs condensate) ip c . The bar 
on the fermion propagator G (x, y) means that it is taken full, while the non-bar notation indicates free propagators, 

as it is the case for the gauge propagator D lil/ (x — y) — J ( ^ ? 4 e '%lJ (sy«/ ~ (1 — «£-ie ) ' anc ^ *he sca l ar one 

A(x — y) = J (27r) 4 E ^+M i -if^ w ^ tn ^ 2 ~ i*^c + M 2 - I n general T 2 (G, y c ) represents the sum of two and higher 
loop two-particle irreducible vacuum diagrams. In the current approximation, as all propagators but the fermion's 
are taken free, 1^ is two-particle irreducible with respect to fermion lines only pq |. In the present weakly coupled 
theory one can use the lowest (two-loop) approximation for T 2 . This corresponds to the so called quenched ladder 
approximation, on which all vertices are taken bare. In this case T 2 is 

r 2 (G,^ c ) = 9 - j 'd 4 xd i ytr[G(x 1 y)rG(y 7 x) 1 »D^(x,y)] 

-Lf d^ytr ( 7 ^G (x, x)) D^(x - y)tr ( 7 "G (y, y)) 

+ ^.J d^xdStr [G (x, y) G (y, x) A(x, y)] 

_^ j d A xd A ytr (G (x, x)) A(x ~ y)tr (G (y, y)) (9) 

The extremum equations (|E|) and (|7|) correspond, respectively, to the Schwinger-Dyson (SD) equation for the fermion 
self-energy operator £ (gap equation), and to the usual minimum equation for the expectation value of the scalar 
field, which in the presence of the magnetic field has to be determined in a self-consistent way, that is, simultaneously 
with the gap equation. 

Although we have introduced a bare scalar mass p, in (0) , because we are interested in the possibility of a dynamically 
generated scalar mass, we take the limit p — > at the end of our calculations. 

The second to fifth terms in the effective action (@) correspond to one-loop contribution. Their evaluation is quite 
straightforward (the scalar self-interaction can be renormalized in the usual way |27|]), with the exception of the 
fermion contributions, which contain the background magnetic field. Then, let us calculate explicitly the one-loop 
fermion contribution coming from the term 

Tf ] =-iTr[G- 1 (i Pc )G\ (10) 

in (^j). Here G (x, y) denotes the full fermion propagator, which can be written as [jl3| 

W dj^K (— i—) E p (y) (11) 

and G~ l (ip c ) is the free fermion inverse propagator in the presence of a constant magnetic field B along the third 
axis, 

G- 1 (x, y, <Pc)=J2[ dP ^ 2 ^ E p (x) (j.p + mo) E p (y) , (12) 



with p = (pojO, —y/2gBk,p3), and mo = \ y <p Cl the fermion mass appearing after the shift ip — > p + <p c in the Higgs 
field. 

In the above equations we have introduced the Ritus' E p functions |]28[ . These orthonormal function-matrices 
provide an alternative method to the Schwinger's approach to problems of QFT on electromagnetic backgroundsn. 



1 For an application of Ritus' method to the QED Schwinger-Dyson equation in a magnetic field see ref. (la]. 



The E p representation is obtained forming the eigenfunction-matrices of the fermion mass operator 

E p (x)=Y,E pa (x)A(a), (13) 

(7 

where 

A(a) =diag(6 rT i,6< 7 -i,d< 7 i,5 cr -i), a = ±1, (14) 

and the E pa functions are given by 

E pa {x) = N(n)e l ^ x ° + P 2x2+ P^D n (p) (15) 

with D n (p) being the parabolic cylinder functions |E9[ of argument p = \[2gB(x\ — -^g-) and positive integer index 

n = n (k,<j) = k+~-~ n = 0,1,2,..., (16) 

and N(n) = (AirgB)^ jyn\ being a normalization factor. Here p represents the set (po,P2,P3,k), which determines 
the eigenvalue p 2 — — pp + p\ + 2gBk in (7^ (id p ~ gA^) 2 i\} p — p 2 ip p (for details and notation see |n| and |l6|]). In 
Eq. ( |T5| ) we are considering the case of a purely magnetic field background (crossed field case) directed along the 
z-direction (without loss of generality we assume that sign(gB) = 1). 
One can easily check that the E p functions are orthonormal 



^ d i xE p ,(x)E p (x) = (27r) 4 5^(p-p') = (2ir) 4 5 kk ,S( Po - p' )5(p 2 - p' 2 )6(p 3 - p' 3 ) (17) 

and complete 

J2 J d 3 P E p (x)E p (y) =J2 I ' d Po dp 2 dp 3 E p (x)E p (y) = (2tt) 4 S^(x - y) (18) 

k k 

Here we have used E p (x) — 7°E , t"'° 



Using Eqs. (O) and (fllh in (M), the last one can be expressed as 

lf> =-%( d 4 xd 4 y E / 7^4 E / P^ Tr i E v (*) (7-P + mo) ^ (2/) 
•/ fc J (27T) k , J (Zir) 

xE ^ y) {wTW))^ ix)} (19) 

Making use of the property (FTTH one can easily integrate in y and p' to obtain 

r?> = -^ ? /|^{ E „ W (^) B , W } « 

At this point we need to consider the structure of the mass operator £ introduced in ref. |L7[ 

E(p) = z » (p)i ■ p„ + z ±(vh ■ p, + m (p) ( 21 ) 



where p = (po,0,0,P3), p ± — (0,0, — \/2gBk 1 Q). The coefficients Z..(p), Zj_(p)and m(p) are functions of p . Notice 
the usual separation in the presence of a magnetic field between parallel and perpendicular variables. Then, taking 
into account (|2l|), the contribution ( P0|) can be written as 

The integral in x yields 



if = -i (2.) 4 S^ (0) V / J^Tr { — — 1 -P^\ TT 1 (23) 

f k-J ( 27r ) \(l + Z ll ) 1 .p+(l + Z 1 _) 1 -p ± +m(p)f 

where the notation S^(k) — S(ko)S(k2)S(k3) is understood. After taking the trace, integrating in p2 and doing the 
Wick rotation to Euclidean coordinates, we obtain 

V) = 8tt 5 .B^ J (0) > / dp 4 dp 3 — - (24) 

We are interested in the contribution of r , to the mininium equations (0) and (ffl). In the case of the gap equation, 
such a contribution can be found directly from Eq. (|10|), differentiating with respect to G. On the other hand, the 
contribution of I, to the scalar vev extremum equation takes the form 

^p(!) p i— \ 

— -L- = 8ir\ v gB6^{0) V / dp 4 dp 3 - 7 ^4= — (25) 

d<Pc v k-J Vl + ( l + z ±) 2 P ± + m2 (p) 

where the solution)^ Z = of the SD equation (g) was explicitly used. 

At large magnetic field, the main contribution to the sum in k comes from the LLL, i.e. p 2 = 2gBk = 0. Then, in 
this approximation Eq. (Eq) becomes 

— J- = 8nX y gBS^(0) / d Pi dp 3 _ 2 ^ _ (26) 

dip c J p|+m 2 (p||) 

In general, the dynamical mass rn(p«) depends on the momentum. However, we can expect that, similarly to QED 



L2[, |18| , m{jpn) behaves as a constant in the infrared region, and diminishes with increasing 



P 



Therefore, the main 



contribution to the integral in Eq. (Eq) will come from the infrared region 



P 



< \/gB. From the above discussion, it 



is reasonable to approximate the function m{jp\\) by a constant solution m{p\\) ~ m(o) — m and use y ' gB as a natural 
cutoff. This leads us to the final result 



.(i) 



dip c 



V^\gBm\n ( ^) = -V^\ y < W> > 
2-k 1 \m ) 



(27) 



where V^ 4 ' represents an infinite four dimensional volume, and < ipip >= iTr {G(x,x)j = — g 27r ? " In I ^ ) denotes 

the fermion-antifermion condensate |13| |, |14| induced by the external magnetic field. 

It is worth to notice the following. It is a well known fact that in the present model, in the absence of a magnetic 

field, the effective action (potential) has a non-zero minimum at the one-loop level, but this minimum lies far outside 

the expected range of validity of the one-loop approximation, even for arbitrarily small coupling constant, so it must 

. dr (1) 

be rejected as an artifact of the used approximation pfj[ . When a magnetic field is present, the term » f , being 

proportional to the fermion condensate, dominates the radiative corrections in the minimum equation (R). Due to 
this, the dynamically generated fermion condensate gives rise to a non zero scalar minimum that is in agreement 
with the used approximation. In other words, thanks to the magnetic field, a consistent minimum solution can be 
generated by radiative corrections. In this sense, a sort of non-perturbative Coleman- Weinberg mechanism takes 
place, with the difference that in the present case no dimensional transmutation is needed. Since the theory already 
contains a dimensional parameter, the magnetic field B, there is no need to include scalar-gauge interactions in order 
to trade a dimensionless coupling for the dimensional parameter ip c f30f. No constraint between the couplings has to 
be assumed, except that they are all sufficiently weak. 



The demostration that Z = is a solution of the gap equation in the present theory can be done along the same line of 
reasoning followed in the Appendix of the first paper of ref. |fq|. 



The two-loop contributions are a little more involved. As we have not enough space in a letter to give all the 
detailed calculations, we will explicitly show, for the sake of understanding, the evaluation of one term. The others 
can be found in a similar way. The complete calculation will be published elsewhere. 

First, notice that the second and fourth term in Eq. (||) generate tadpole diagrams in the SD equation (g). It is 
easy to realize that the tadpole diagram with the gauge-fermion vertex vanishes. However, the tadpole associated to 
the scalar-fermion vertex is not zero when ip c ^ and, as shown below, it has a significant contribution to the gap 
equation. Let us evaluate this tadpole contribution, which we denote by ^ . 



E T (^ y) = i 5 -^ = -t\y(x -y) J d**A(a: - z)tr [g(z, z) 
We can transform Eq. (|23) to momentum space with the help of the E p (x) functions to obtain 



J d i xd i yE p {x)Y. T {x,y)E pl (y) = (2i:f^(p - p') £ T (p) 



(28) 



— , d 4 n P ig( x - Z ) 

-i\ 2 y / d A xd A zE p (x) l ' 



*£/f^M 



(2tt) 4 q 2 + M 2 + i e 
1 



— E p « (z) E p , (x) 



Taking into account that [[L3| 

J d 4 xe lqx E p {x)E p , (x) = (2irf6 (3 \p' + q-p)e tqi ^ +P2 ^ 29B e-^ /2 

e i(n-n')<p 

one can integrate in x and z to find 

(2^)^\p-p')Y J T {p) = ~^l f AE [d 3 p"5^(g)S^{p' + q -p) 



(29) 



(30) 



xe 



-2 e *9i(P2+P2-2p 2 ")/2g-B ^^ e i(n-n')<p 



q 2 + M 2 + ie ^ y/n(k,<T)\n'{k',a)\ 



J„„'(§1)A(<t) 



£ \Tr 



A K)^T 



7-p" + E(p") 



n"(fe",a")! Jn """ (5l)r ' 



(31) 



This equation can be further simplified after taking the trace and using the small ^approximation of the J— functions 



[max(n,n")]! |„_ w »| 



(32) 



which can be justified by the presence of the exponential factor e~ 9 ^ in the integrand of Eq. (plj). Moreover, thanks 
to the delta 5^ 3 '(q), the integrations in qo, q2, qs are trivial. Thus, from the previous considerations and using the 
properties of the A matrices [ [13[ , we arrive at 

. . _t r , , r -~i P i( ?i(P2+P2— 2p2")/ 2 s-B 

(2tt) V 4 ) (p ~ p>) J] m = -2i\ 2 y 8^ (p> -p)j d qi J2J d 3 P " e-o- g 2 + M2+ . e 



2m(p") 



(l + z) p "2 + (l + Z ± ) 2 ^"2 + m2(p") 



(33) 



Finally, after integrating in q± and P2" , and transforming to Euclidean space, we get 
T A 2 gB ^ f „ „ m(p") 

E =2^Jp2^ dp4P3 7 x 3 2 (34) 

V 11/ II 

Similarly to the analysis done in the one-loop case, it can be seen that the main contribution to Eq. (J34J) comes 
from the fc" = term of the sum. Taking into account this and using the solution Z = 0, the tadpole contribution 
to the gap equation reduces to 

A 2 nB f m (P ") 

^ = 2^JP J dpldp " 3 p»* + m*(p») (35) 

11 11 

One can recognize here the same integral that led us to the fermion condensate in Eq. (|2^). Therefore, we find 

Note that the tadpole term is proportional to the magnetic field and inversely proportional to the scalar mass. This 
functional dependence will be responsible for the notorious increment of the mass solution in this model as compared 
to other theories, as discussed below. 

Taking into account the leading contributions to Eqs. (g) and (0) at large magnetic field, one arrives at the following 
minimum equations for the fermion mass and the scalar vev respectively, 

— (^!)£'» 2 (S)+iJr-KS) <-) 

Eq. (Bq) can be further simplified by noting that at large B one can neglect the terms ~ A 2 coming from the 
one-loop scalar self-interaction, compared to the term coming from the fermion condensate contribution ~ < ifnp >. 
Then, the scalar minimum satisfies 

,„3 _, \ 3; 

A 7r z V m z 



Vl^¥-£rnln( 9 4) (39) 



Similarly, the gauge and scalar bubble diagram contributions to the gap equation, (second term in (|37|)), are 
igligible compared to the tadpole contribution, (third term in (|37|)), so the equation becomes 



m~mo + -^-^gBm\n( ^ ) (40) 



Substituting @ in Eq. @ , it is found 



where the coefficient k satisfies 



m~-^=y/g~B (41) 



k1hk~1.4— (42) 

\y 



The corresponding solution for the Higgs vev is 



^ik,^~ B < 43 > 



It is easy to check that for weak couplings the solutions ( |4l| ) and ([43]) are in agreement with the used approximations. 
Note that there is no zero solution for the scalar vev in this large field approximation. Both, the minimum of the scalar 
field and the dynamically generated mass are driven by the external magnetic field. As the dynamical fermion mass 
breaks the discrete chiral symmetry (|5j) , this model might be considered as one more example of the phenomenon 
of magnetic catalysis. However, if the current model is extended to include complex scalars (complex scalars do not 
change at all the conclusions of this paper) and non-simple gauge groups, as 5/7(2) x U(l) in the electroweak model, 
the symmetry breaking phenomenon will have a different nature. There, we do not have chiral symmetry, but the 
magnetic- field-driven scalar minimum will break the gauge symmetry by giving mass to the gauge fields coupled to it. 
Thus, as we claimed at the beginning of the paper, in richer models with scalar fields, the magnetic field can either 
reinforce gauge symmetry breaking if it is already broken, or induce it, through non-perturbative effects. 

Comparing the induced fermion dynamical mass Eq. (uW with the mass generated when no scalar field is present 
Ufl , fll3|| , m ~ \/\gB\ exp [— 2n/g] , or with the one obtained when the couplings are fine-tuned so the scalar vev is set 



to zero Q , m w \J\gB\ exp 



-2ttW1/ [g 



,2 



one realizes that the scalar interactions, when taken into account 



in a self-consistent way, can dramatically strengthen the generation of mass. This observation is easy to corroborate 
by direct evaluation of the mass (Ell) for typical values of Yukawa coupling X y and scalar self-coupling A. For instance, 
if we take \ y = 0.7, which is the approximate value of the Yukawa coupling for the top quark, and A = 0.4, which 
corresponds to a Higgs mass of 115 Gev, we find m ~ 0.6y|^B|. The same Yukawa coupling, on the other hand, gives 
just m ~ 10~ 5 \Z\g~B\ if the scalar vev is fine-tuned to zero |l6f . In the case of pure gauge- fermion interactions, as for 
instance in QED, the generated mass would be even much smaller p2| , jig] . 

A non-simple group extension of the model discussed in this work could be of interest as an effective theory in 
condensed matter problems, where SU(2) x U(l) gauge theories (although without Higgs fields) have been previously 
proposed to describe the rich phase structure of high T c superconductors |21j] . 

It seems, however, that the most immediate physical extension of the present model would be the electroweak 
theory. This case is particularly interesting in the light of the recent works on the role of magnetic fields in electroweak 
baryogenesis M- finl. In view of the magnetic-field-driven non-perturbative enlargement of the Higgs vev, Eq. (Ejj), 
it is possible that in the electroweak theory the increase in the scalar vev will be large enough to guarantee the 
baryogenesis condition (pi) . It remains therefore as an open question whether the effect found in this paper can 
influence the recent conclusions ||- |f(|] about baryogenesis in the presence of primordial magnetic fields. 
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